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Abstract. We prove that the q-states Potts model on graph is spontaneously 
magnetized at finite temperature if and only if the graph presents percolation on the 
average. Percolation on the average is a combinatorial problem defined by averaging 
over all the sites of the graph the probability of belonging to a cluster of a given size. 
In the paper we obtain an inequality between this average probability and the average 
magnetization, which is a typical extensive function describing the thermodynamic 
behaviour of the model. 
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1. Introduction 

The interplay between spin models and percolation has been put into evidence since the 
fundamental work of Fortuin and Kasteleyn |L, where it is shown that percolation is 
the limit for q — > 1 of the random cluster representation of q-states Potts model. This 
representation, as already pointed out in [1J, is very general and can be performed on 
any discrete structure, i.e. on graph. 

A further step has been the result proven in |2j, showing that, on lattices, the 
existence of percolation implies Potts transition for any value of q. The main purpose 
of this work is the extension of the result to generic networks. 

An important issue in this direction has been a mathematical paper [3] proving 
that Potts model has more than one Gibbs measure if it is defined on graphs showing 
percolation. The multiplicity of Gibbs measures is the most used definition of symmetry 
breaking in mathematical literature. However, it does not always correspond to the 
thermodynamic concept of phase transition. 

Indeed, in the study of the thermodynamic properties of physical systems, we are 
interested in the free energy and its derivatives such as the average magnetization, the 
susceptibility, the specific heat, which are typical average extensive quantities. On the 
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Figure 1. The brush graph is obtained by adding a linear chain to each site of a 
two-dimensional lattice. 

other hand, the existence of a single Gibbs measure, as explained in [3], depends on a 
local parameter, that is the magnetization of a site of the graph. On lattices, translation 
invariance implies that local and average quantities coincide (the proof of |2j depends 
on this fact). However, for inhomogeneous structures the results can be different in the 
two cases. For example, on the brush graph (see Figure 1) as we will prove in detail, 
there is not a single Gibbs measure, but the thermodynamic functions do not present 
discontinuities and the spontaneous magnetization is zero. 

A generalization of the result [2j, which takes into account the thermodynamic 
behaviour of the system, requires the introduction of the new problem of percolation 
on the average. A graph G is said to present percolation on the average if the average, 
over all sites, of the probability of belonging to a cluster of size greater than / remains 
positive when I — > oo. We will prove that the q-states Potts model on a graph G is 
spontaneously magnetized if and only if G presents percolation on the average. 

Therefore, the study of spin models with discrete symmetry leads to a classification 
of graph topology which can be formulated in terms of a purely combinatorial problem. 

An analogous result 4] has been previously obtained for classical spin models 
with continuous symmetry, where the thermodynamic behaviour depends on a different 
combinatorial problem, i.e. random walks. In particular, the models present 
spontaneous magnetization if and only if the graph is transient on the average jH] (i.e. a 
graph where the average probability for a random-walk of ever returning to the starting 
site is smaller than 1). 

This two classifications are different, indeed, it is well known that lattices presenting 
percolation (i.e. lattices with dimension d > 2) are transient on the average only if d > 3. 
In recent works t3j, it has been proved that the "anomalous" behaviour of presenting 
average percolation in a recurrent on the average structure is typical not only of two 
dimensional lattices, but also of low dimensional networks, such as the Sierpinski carpet. 
Therefore, a general result clarifying the relation between the two classifications is still 
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lacking. Furthermore, up to now we do not know if percolation on the average is related 
to graph topology by a simple parameter, as it is the spectral dimension |7| in the case 
of transience on the average. 

The paper is organized as follows. In the next section we introduce the Potts model 
and percolation on graph, we define percolation on the average and we state the main 
result of the paper, i.e. Potts model is spontaneously magnetized if and only it is defined 
on a graph presenting percolation on the average. In the last 2 sections we prove the 
theorem. 



2. The q-states Potts model and percolation on graph 

A graph G (see e.g. |S] for an introduction to graph theory) is a countable set 
V of vertices or sites i connected pairwise by a set E of unoriented edges or links 
= Connected sites are called nearest neighbours and we denote with Zi the 

connectivity (number of neighbours) of the site i. A path is a sequence of consecutive 
links {(i, k)(k, h) . . . (n,m)(m, j)}. The chemical distance is the length (number of 
links) of the shortest path connecting the sites i and j. The Van Hove sphere S or of 
center o and radius r is the set of the sites of G such that r 0) j < r. We will call E 0>r the 
set of all E E such that i,j E S 0>r , dS r the set of the vertices of S r connected with 
the rest of the graph and we denote with | • | the cardinality of a set. Let (pi be a real 
function of the sites of an infinite graph, the average in the thermodynamic limit is 

E & 

0= lim^— . (1) 

r— >oo ,S 

I '-'OjT | 

Interesting properties of the thermodynamic average, such as the independence of the 
choice of the center o, are proved in j^J. Therefore, in the following we will drop the index 
o. The measure \\S\ \ of a subset S of V is the average value x(S) °f ^s characteristic 
function Xi(S) defined by Xi(S) — i H i E S and Xi{S) = if i E~ S. 

Important constraints on graph topology follow from the requirement to describe 
a physical network. Real systems, indeed, have bounded local energy and they are 
embedded in a finite dimensional space. Therefore, we will consider only connected 
graphs with bounded connectivity (zi < z max Vz E E) and such that: 

lim M = (2) 

r -*°° \S r \ 

The definition of the q-states Potts model (see e.g. for a review) on an infinite 
graph G first requires to introduce the model on the Van Hove spheres S r . For each site 
i E S r Si is a q state function Sj = 1 . . . q and the Hamiltonian is given by: 

H r = ]T (l-6( Si , Sj )) + hJ2(l-S(si,l)) (3) 

{i,j)eE r i£S r 

where the Kronecker delta-function is 5(si, Sj) = 1 if Sj = Sj and 6(si, Sj) = otherwise. 
Notice that the case q = 2 is the Ising model. In the canonical ensemble the partition 
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function Z r is given by the sum of the Boltzmann weight exp(—j3H r ) over all possible 
configurations {sj}; (3 represents the inverse temperature of the system. Thermodynamic 
properties of statistical models are described by extensive order parameters, in this case 
the average magnetization: 

M r (P, q, h) = \Sr\- 1 (S(si, 1) - q' l )r = 

\Sr\- X T, Z r l E Ws t , 1) - q-^e-M . (4) 

ies r S 1 ...S Sr 

where (-) r denotes the thermal average. The Potts model on G presents spontaneous 
magnetization if it exists (3 C such that for all (3 > (3 C 

lim M(3, q, h) = lim lim MJ/3, q, h) > . (5) 

The existence of the thermodynamic limit r — > oo is always assumed in this article. 

Percolation (see e.g. [TOi m] for a mathematical and a physical treatment) is defined 
by introducing a probability p, < p < 1 and each link G E is declared to be 

open with probability p and closed with probability 1 — p independently. The cluster 
Ci containing the site i is the set of all the sites connected with % by a path of open 
links. The i-size of the cluster Cj is the maximum chemical distance from % of the sites 
of Ci. We call Pi(l,p) the probability for i to belong to a cluster of i-size > I. A 
graph G is said to present (local) percolation if it exists a probability p < 1 such that 
lim^oo Pi(l,p) > 0. On a graph with bounded connectivity this property is independent 
of i [IDj. 

A graph presents percolation on the average if it exists a probability p such that 

lim P(l,p) = lim lim |5 r | _1 V Pi(l,p) > . (6) 

tes r 

Percolation on the average is a new combinatorial problem and it gives rise to a graph 
classification which turns out to be fundamental for understanding the thermodynamic 
behaviour of physical models on graph. In the main theorem of the paper, indeed, we 
prove that the q-states Potts model is spontaneously magnetized if and only if it is 
defined on a graph presenting percolation on the average, obtaining a generalization to 
inhomogeneous structures of the classical result for lattices by Aizenman et al. .2]. The 
proof is mainly inspired by [3] and jT] . 

First, following pQ, we introduce a representation of the q-states Potts model in 
terms of the q-random cluster model on the supplemented graph G'. Then, from an 
important property of the random cluster stated in t 3], we get an inequality between 
the local magnetization and the percolation probability on the supplemented graph: this 
is an extension of [3] to the case h ^ 0. In the second part of the proof we take the 
thermodynamic average and we show that the physical requirements on the graph allow 
to formulate the inequality for the average magnetization in terms of average percolation 
on G. The theorem directly follows from this inequality. 

Let us first show that local and average percolation define different classifications 
of inhomogeneous networks. In particular, the brush graph (Figure 1) presents local 
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percolation but no percolation on the average. If we choose p larger than the threshold 
of the two-dimensional percolation and we call r the distance of % from the plane 
of the brush, we have that Pi(l,p) > P sq (l — r,p)p r (P sq (l,p) is the probability 
that a site of a two-dimensional lattice belongs to a cluster of size > I) and then 
lirn^oo Pi(l, p) > p r lirn^oo P sq (l — r,p) > 0. On the other hand for all the sites of 
the brush at a distance from the plane greater than I we have Pi(l,p) = 2p l — p 21 , 
2p l — p 21 is the probability for a site of a linear chain to belong to a cluster of i-size 
greater than /. We call Ri the subset of the brush graph given by the sites whose 
distance from the plane is smaller than /. We have =0 and 

lim P(l,p) < lim ||^|| + (1 - \ \R l \\)(2p l - p 21 ) = . (7) 

Therefore, on the brush graph the q-states Potts model does not have a single 
Gibbs measure however the system is not spontaneously magnetized and the 
thermodynamic functions are analytical. 

3. Random cluster representation of the Potts model 

Let us define the supplemented q-random cluster model [I]. We call G' the graph 
obtained by adding to G a supplementary site o connected by a new link to each site of 
G. The supplemented spheres are obtained by adding to the spheres S r the new vertex o 
and the relevant edges. We will call S' r and E' r the set of sites and the set of links of the 
supplemented sphere. The configurations £ of the supplemented random cluster model 
on S' r are obtained by declaring each link of E' r open or closed. The supplemented 
q-random cluster model pQ is defined by choosing each configuration according the 
probability distribution: 

Mp,Po,9,0 =P Ns{0 ■ (1 -p)W~ N *® (0 • (1 - Po frh»o(V q c(0 (8) 

where p, p a and q are real parameters of the model (<7>l,0<p<l and < p Q < 1. 
Ng(£) is the number of open links in E r , N (£) the number of open links connecting o to 
S r and C(£) the number of clusters in the configuration (|5V| and \E r \ are respectively 
the cardinality of S r and E r ). For q = 1 and p Q = the random cluster model is exactly 
the percolation model restricted to the sphere S r . 

As already pointed out in the original work by Fortuin and Kasteleyn lj, q-states 
Potts model and q-random cluster model are equivalent on any graph G. Indeed if we 
fix p = 1 — e _/3 and p = 1 — e~ /3fc , their partition functions coincide and 

(5(s t ,l)-q- 1 ) r = P t °(p,p ,q)^ J F?{S)»(p, Po ,q,S) (9) 

where F°(£) — 1 if o G Cj and F°(£) = otherwise. Then, P°(p,p ,q) represents the 
probability for the site i to belong to the same cluster of the supplemented site o. 

The space of configurations of the random-cluster model is equipped by the partial 
order: £ ^ 77 if the set of open links in £ is included or equal to the set of open links 
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in \x. We will say that the probability distribution stochastically dominates the 
probability distribution //(£) {p! -<d p) if ; 

/ fdy! < J fdfi (10) 

for all increasing function /. 

In |3j it is proved that important inequalities for the probability distributions follow 
from FKG J2| theorem. These inequalities can be easily generalized to the case of 
supplemented graph, obtaining: 

Kp'iP'oAiO ±d fJ>(p,Po,q,£) —D /x(p,p o ,l,0 (11) 

where p' = p(q — p(q — l))^ 1 and p' Q = p Q (q — p (q — l)) -1 . 

Since F°(£) is an increasing function, from (jHJ) and (jllj) we have: 

P°(p',p' , 1) < (5(s h 1) - q- l ) r < P t °(p, Po , 1) (12) 

(fT2"j) is an extension to the case h ^ of the inequality for local magnetization proved 



4. Spontaneous magnetization and percolation on the average 

Let us first take the average of ((E)) over all the sites of the sphere S r . We obtain 

\Sr\- 1 mp',P' A) < M r ((3,q,h) < \Sr\- 1 2 /f(p,p ,l) (13) 

Inequalities (|T3*j) provide upper and lower bounds for the magnetization of the Potts 
model in terms of the percolation model defined on the supplemented Sphere S' r . Let 
us reformulate these bounds in terms of the percolation on G. 

First we prove a property of graphs with bounded connectivity and satisfying (j2J. 
In the probability measure defined on S r by (JHJ with q = 1 and p a = 0, we call P itr (l,p) 
the probability for the site i to belong to a cluster of i-size greater than I. Let us show 
that in the thermodynamic averages one can substitute P{ jr (l,p) with the probability 
defined on the infinite graph Pj(Z,p). We will call Sf r the subset of S r of all the sites i 
such that the distance of i from the border dS r is larger than I and Sf r its complementary 
in S r . For the sites i G Sf r all the clusters Cj of i-size smaller than / are included in S r , 
then the probability of belonging (or not belonging) to one of these clusters is equal in 
percolation on G and in percolation on S r . Furthermore we have that: 

||^- || = hm \S r \- l \Si; r \ < hm l^r 1 ^^ = (14) 

where we used (J2J and the boundedness of the coordination number. Then P ijr (l,p) and 
Pi(l,p) differ only on a set of zero measure and in the thermodynamic averages we can 
exchange one for the other. 

Let us now prove the main theorem of the paper. From the boundedness of the 
correlation number and the independence of percolation probabilities, one gets the 
inequality: 

P°{ P ,Po, 1) < (l - (1 - Pof"-) (1 - Pi,r(l>P)) + P U^P) ■ (15) 
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Indeed (z l max — l)/(z max + 1) is the maximum number of sites in a cluster of i-size < / 
and then (1 — (1 — p ) zImax ) is an upper bound on the probability that a site of the cluster 
is connected to o. From (jlHjl and (fTK|) taking the thermodynamic limit r->oowe have: 

M(P,q,h) < lim (f 1 " (1 ~ VoY La * ) (1 - Pi{l,p)) + W,P)) 

Letting h — > 0, since also p a — > 0, we get 

lim M(p, q, h) < lim ( 1 - (1 - p Y™A + P(l,p) = P(l,p) . (17) 

Now we have to consider the first inequality in f)13|) . From the probability independence 
in percolation, we get 

P°(p', P ' , 1) > (l - (1 P*,r(i,l/) • (18) 

Here / + 1 is the minimum number of sites in a cluster of i-size > /, and then 
(1 — (1 — p' ) l+1 ) is a lower bound on the probability that a site of the cluster is connected 
to o. Taking the thermodynamic average we have 

M(/3, q, h) > (l - (1 - p' ) l+1 ) ■ (19) 

Letting first / — > oo and then /i^Owe have 

lim M(P, q, h) > lim P(l,p') (20) 
Finally, from (|17p and (|2U|) one has 



lim P(l, p') < lim M{(3, q, h) < lim P(l, p) (21) 

l— too h^O l^oo 

If the graph presents percolation on the average, it exists a probability p' such that 
lim^oo P(l,p') > and then the first inequality in (|2*T|) implies that the model for 
(3 = log(l +p'q(l —p')^ 1 ) is magnetized. On the other hand, if lim^oo P(l,p) = for all 
values of p, the second inequality in (|21j) implies that the magnetization of the system 
is zero for all temperatures and this concludes the proof. 

The result can be easily generalized to Potts models with disordered ferromagnetic 
couplings described by the Hamiltonian H r = J2(i,j)eE r JijO- ~ $(. s ii s j)) + hJ2i^s r (^ ~ 
5(si, 1)) with < e < Jij < K for all (ij) E G. 

On the other hand an important open question is the possibility to extend the proof 
to different ferromagnetic spin models with discrete symmetry, such as the clock model. 
As mentioned in the introduction, another basic problem is a better understanding of 
the relation between percolation on the average, transience on the average and graph 
topology. 
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